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ABSTRACT 

One of the most important and poorly-understood issues in structure formation is the role of outflows 
driven by active galactic nuclei (AGN). Using large-scale cosmological simulations, we compute the 
impact of such outflows on the small-scale distribution of the cosmic microwave background (CMB). 
Like gravitationally-heated structures, AGN outflows induce CMB distortions both through thermal 
motions and peculiar velocities, by processes known as the thermal and kinetic Sunyaev-Zel'dovich 
(SZ) effects, respectively. For AGN outflows the thermal SZ effect is dominant, doubling the angular 
power spectrum on arcminute scales. But the most distinct imprint of AGN feedback is a substantial 
increase in the thermal SZ distortions around elliptical galaxies, post-starburst ellipticals, and quasars, 
which is linearly proportional to the outflow energy. While point source subtraction is difficult for 
quasars, we show that by appropriately stacking microwave measurements around early type galaxies, 
the new generation of small-scale microwave telescopes will be able to directly measure AGN feedback 
at the level important for current theoretical models. 

Subject headings: cosmic microwave background - galaxies: evolution - quasars: general - intergalactic 
medium - large-scale structure of the universe 



1. INTRODUCTION 

Large-scale measurements of cosmic microwave back- 
ground (CMB) anisotropies have played a central role in 
the development of the modern cosmological model. By 
giving us a picture of the universe during the epoch of 
linear fluctuations, they have provided us with a foun- 
dation from which to interpret the subsequent growth 
of structure. Thus from the spectacular initial detec- 
tions made by the COBE satellite (Smoot et al. 1992), 
to the detailed flatness constraints provided by balloon- 
born missions (Hanany et al. 2000; Mauskopf et al. 2000), 
to the precision measurements of cosmological parame- 
ters derived from the WMAP observations (Spergel et 
al. 2003; 2007), our understanding of structure forma- 
tion has progressed hand-in-hand with improvements in 
large-scale CMB anisotropy measurements. 

On angular scales below w 10', the science potential 
of CMB measurements has yet to be realized. On these 
scales, Silk damping washes out the primary anisotropies, 
while a host of smaller-scale secondary anisotropies are 
imprinted as the photons propagate across the universe. 
One of the most important sources of these effects is the 
scattering of CMB photons by hot electrons, a process 
that was first studied by Sunyaev & Zel'dovich (1970; 
1972). This process can be divided into two contri- 
butions. The largest of these is the thermal Sunyaev- 
Zel'dovich (SZ) effect in which inverse Compton scat- 
tering preferentially increases the energy of CMB pho- 
tons passing through hot and dense regions, introducing 
anisotropies with a distinctive frequency dependence. A 
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smaller contribution arises from the kinetic SZ effect in 
which Doppler scattering in dense regions with signifi- 
cant peculiar motions induces fluctuations with the same 
frequency dependence as the CMB itself. 

The most important source of such hot and dense re- 
gions is undoubtedly the gravitational collapse of gas into 
large dark matter halos, and numerous numerical simu- 
lations of the Sunyaev-Zel'dovich effect from this process 
have been carried out by several groups (e.g. Scaramella, 
Cen, & Ostriker 1993; Hobson & Magueijo 1996; da Dilva 
et al. 2000; Refregier et al. 2000; Seljak, Burwell, & Pen 
2001; Springel, White, k. Hernquist 2001; Zhang, Pen, 
& Wang 2002; Roncarelli et al. 2007). Yet dark-matter 
driven gravitational heating need not provide the domi- 
nant SZ contribution at all scales and in all environments. 
In fact the gas responsible for the largest SZ distor- 
tions, the intracluster medium (ICM) in galaxy clusters, 
is observed to have been substantially heated by non- 
gravitational sources (Cavaliere et al. 1998; Kravtsov & 
Yepes 2000; Wu et al. 2000; Babul et al. 2002). In galaxy 
groups, such nongravitational effects are even more se- 
vere, causing groups with a large range of X-ray lumi- 
nosities to all be heated to a similar value of « 1 keV per 
baryon (e.g. Arnaud & Evrard 1999; Helsdon & Ponman 
2000). 

But perhaps the most dramatic, yet poorly understood 
manifestation of nongravitational heating is in establish- 
ing the observed "downsizing" (Cowie et al. 1996) trend 
in the evolution of star-forming galaxies and active galac- 
tic nuclei (AGN). Here the issue is that since z rj 2 
the characteristic mass scale of star-forming galaxies and 
the typical luminosities of AGN have dropped by over 
an order of magnitude (e.g. Arnouts et al. 2005, Treu et 
aZ.2005; Pei 1995; Ueda et al. 2003; Barger et aZ.2005). 
This "antiheirarchical" trend is in direct conflict with the 
expectations of the long-standing model of structure for- 
mation in which gas condensation and heating is driven 
purely by dark mater halos, which grow hierarchically 
by accretion and merging over time. In this picture, as 
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gas falls into potential wells, it is shock heated and must 
radiate this energy away before forming stars (Rees & 
Ostriker 1977; Silk 1977). The larger the structure, the 
longer it takes to cool, and thus the formation of history 
of galaxies is even more hierarchical than the dark matter 
history. 

Although early modeling did not reproduce this down- 
sizing behavior, the trend has been successfully repro- 
duced in more recent theoretical models that include a 
large source of heating associated with outflows from 
AGN (Scannapieco & Oh 2004; Binney 2004; Granato 
et al. 2004; Scannapieco, Silk, & Bouwens 2005; Di Mat- 
teo, Springel, & Hernquist 2005; Croton et al. 2006; Cat- 
taneo et al. 2006; Thacker, Scannapieco, & Couchman 
2006, hereafter TSC06; Di Matteo et al. 2007). In this 
picture, AGN outflows associated with broad-absorption 
line winds and radio jets heat the surrounding intergalac- 
tic medium (IGM) to sufficiently high temperatures to 
prevent it from cooling and thus from forming further 
generations of stars and AGN. This feedback requires an 
energetic outflow, driven by a large AGN, to be effective 
in the dense, high-redshift IGM, while in the more ten- 
uous low-redshift IGM, equivalently long cooling times 
can be achieved by less energetic winds. The lower the 
rcdshift, the smaller the galaxy that is able to exert effi- 
cient feedback, thus resulting in cosmic downsizing. 

However, the details of AGN feedback remain ex- 
tremely uncertain. Heating may be impulsive (e.g. 
TSC06), more gradual (e.g. Briiggen, Ruszkowski, & 
Hallman 2005; Croton et al. 2006), or modulated primar- 
ily by the properties of the surrounding material (e.g. 
Dekel & Birnboim 2006). Furthermore, direct measure- 
ments of the kinetic energy input from AGN are notori- 
ously difficult and range from ~ 1% or less (de Kool 2001) 
to ~ 60% of the AGN's total bolometric energy (Char- 
tas et al. 2007) . Finally, several alternative ideas have 
been suggested to explain downsizing outside of the con- 
text of AGN feedback altogether (e.g. Keres, et al. 2005; 
Khochhfar & Ostriker 2007; Birnboim, Dekel, & Neistein 
2007). 

Here we show that the definitive measurement resolv- 
ing this issue may again come from the microwave back- 
ground, through observations of SZ distortions. As it 
directly probes the spatial distribution of heated gas, 
SZ detections are able to place constraints on the pri- 
mary theoretical uncertainty in AGN feedback models, 
namely the nature and degree of nongravitational gas 
heating. Fortunately, the rise in importance of this is- 
sue for theory is paralleled by recent advances in ob- 
servation. Motivated primarily by deriving cosmologi- 
cal constraints though the detection of galaxy clusters 
(e.g. Holder et al. 2000; Majumdar & Mohr 2003; Battye 
& Weller 2003; Schulz & White 2003), a large number 
of blank-field small-scale CMB surveys are now under- 
way, using telescopes that will push into the interest- 
ing regime for AGN feedback. These include surveys 
with the Atacama Cosmology Telescope (Kosowsky et 
al. 2006), the South Pole Telescope (Ruhl et al. 2004),' the 
Atacama Pathfinder Experiment (Dobbs et al. 2006), and 
the Sunyaev-Zel'dovich Array (Loh et al. 2005), many of 
which will be coordinated to overlap with optical sur- 
veys. These advances are particularly important as less 
sensitive small-scale CMB surveys have already hinted 
at an excess of SZ distortions (Mason et al. 2003; Bond 



2005; Dawson et al. 2006; Kuo et al. 2007). 

In this work we make use of the first large-scale cosmo- 
logical hydrodynamical simulation to explicitly include 
AGN feedback (TSC06) to construct detailed maps of 
their imprint on the CMB. Analyzing these maps and 
comparing them with other observables, we are able to 
determine the most promising approaches to using these 
ongoing surveys to constrain AGN feedback. Previous 
studies of nongravitational heating on the SZ background 
focused on the impact of starburst winds at high redshift 
(Majumdar, Nath, & Chiba 2001; White, Hernquist, & 
Springel 2002) and signatures from Population III stars 
(Oh, Cooray,& Kamionkowski 2003). Recently, Chatter- 
jee & Kosowsky (2007) made analytic estimates of the 
impact of quasar winds on the CMB power spectrum and 
the cross-correlation of the SZ effect with optical sources. 
Here we show directly from simulations that while basic 
quantifiers such as the power spectrum are difficult to in- 
terpret, approaches based on the cross-correlation of SZ 
distortions with optical observations will soon provide a 
clean and direct probe of AGN feedback. 

The structure of this work is as follows. In §2 we de- 
scribe our numerical simulations and the methods used to 
construct maps of the CMB distortions from the kinetic 
and thermal SZ effects. In §3 we compute the contri- 
bution of AGN outflows to the CMB power spectrum, 
calculate their impact on the regions surrounding indi- 
vidual galaxies and AGN, and quantify the sensitivity 
of CMB-galaxy cross correlations to AGN feedback. Fi- 
nally, our conclusions are summarized in §4. 

2. METHOD 

2.1. Simulations 

In order to disentangle the SZ signatures of AGN feed- 
back from distortions due purely to gravitational heating, 
we made use of two simulations, the large AGN feed- 
back simulation first presented in TSC06 (using 2 x 640 3 
particles) and a smaller-scale comparison simulation (us- 
ing 2 x 320 3 particles) in which structure formation, gas 
cooling, and star formation were tracked exactly as in 
the fiducial run but no outflows were added. In both 
cases, based on wide range of cosmological constraints 
(e.g. Spergel et al. 2003; Vianna & Liddle 1996; Riess et 
al. 1998; Perlmutter et al. 1999), we adopted Cold Dark 
Matter cosmological model with parameters h = 0.7, flo 
= 0.3, Q A = 0.7, Q b = 0.046, <7 8 = 0.9, and n = 1, where 
h is the Hubble constant in units of 100 km s _1 Mpc -1 , 
i7o, f^A, and Qb are the total matter, vacuum, and bary- 
onic densities in units of the critical density, cr| is the 
variance of linear fluctuations on the 8ft. -1 Mpc scale, and 
n is the "tilt" of the primordial power spectrum. Note 
that these parameters are slightly different than those 
preferred by the more recent WMAP data (Spergel et 
al. 2007) , which was released after our large AGN feed- 
back run had already been completed. Initial conditions 
were computed using the Eisenstein & Hu (1999) transfer 
function. 

In both runs, as in our previous work (e.g. Scanna- 
pieco, Thacker, & Davis 2001), simulations were con- 
ducted with a parallel OpcnMP based implementation of 
the "HYDRA" code (Thacker & Couchman 2006) that 
uses the Adaptive Particle-Particle, Particle-Mesh algo- 
rithm to calculate gravitational forces (Couchman 1991), 
and the smooth particle hydrodynamic (SPH) method to 
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calculate gas forces (Lucy 1977; Gingold & Monaghan 
1977). As the details of this code and our outflow im- 
plementation are described in detail elsewhere (Scanna- 
pieco, Thacker, & Davis 2001; TSC06) here we only sum- 
marize the aspects most relevant to the SZ effect. 

Our study is targeted to relatively large and late form- 
ing structures, and for both runs we have kept the metal- 
licity constant at Z = 0.05, to mimic a moderate level 
of enrichment. Similarly, because the epoch of reioniza- 
tion is poorly known and because reionization has little 
impact on mass scales greater than 10 9 M Q (Barkana, & 
Loeb 1999), we did not include a photoionization back- 
ground in either simulation. In our AGN feedback run, 
which was designed to cover a large range in the AGN 
luminosity function, we used a simulation box of size 
146ft,- 1 Mpc filled with 2 x 640 3 particles, which corre- 
sponded to a dark-matter particle mass of 1.9 x 10 8 M Q 
and a gas particle mass of 2.7 x 10 7 M Q . Our comparison 
simulation, which was carried out purely for the pur- 
poses of the current study, used 320 3 particles in a box 
73/i _1 comoving Mpc on a side, corresponding to the 
same particle masses as in the AGN run. Both runs were 
terminated at z = 1.2, at which point integration was 
becoming expensive due to the single-stepping nature of 
the Hydra code, and impractical in a shared queue en- 
vironment. This means that our results do not contain 
the SZ contribution from lower redshifts, during which 
most galaxy clusters are formed. However, as z = 1.2 is 
well past the peak epoch of AGN activity (eg. Ueda et 
al. 2003; Barger et al. 2005), our results should do well at 
quantifying the AGN contribution to the SZ effect, which 
is our main focus here. 

As in TSC06, bright quasar-phase AGN are associated 
with galaxy mergers, which are tracked by labeling gas 
particles and identifying new groups in which at least 
30% of the accreted mass does not come from a sin- 
gle massive progenitor. Once a merger has been identi- 
fied, we compute the mass of the associated supermassive 
black hole, Mbh, from the circular velocity of the rem- 
nant, v c using the observed Mbh — v c relation (Merrit 
& Ferrarese 2001; Tremaine et al. 2001; Ferrarese 2002) 
which gives 



Here v c is estimated as 
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where G is the gravitational constant, p v (z) is the virial 
density as a function of redshift, and r v is the implied 
virial radius for a group of N gas particles with mass m g 

_ r Nm^noM 1/3 r * 

?v " [ 4/37rp v (z) J ' [6) 

Following Wyithe & Loeb (2002) and Scannapieco & Oh 
(2004), we assume that for each merger the accreting 
black hole shines at its Eddington luminosity (1.2 x 10 38 
ergs s _1 M _1 ) for a time taken to be a fixed fraction 
of the dynamical time of the system, td = 0.055r v /w c = 
5.8 x l^ z n{z)- 1 /' 2 H{z)- 1 . In TSC06 it was demon- 
strated that, apart from a discrepancy for the very lumi- 
nous end, this simple model does extremely well at re- 
producing the observed AGN luminosity function as well 
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Fig. 1. — Evolution of the mass- averaged gas temperature in 
our AGN feedback (solid) and comparison (dashed) runs. 



as the large and small-scale clustering of AGN over the 
full range of simulated redshifts. Furthermore, the dis- 
crepancy at the very luminous end for the simulation as 
compared to the semi-analytic model could be attributed 
to the relative efficiency of shock heating on substructure 
(e.g. see Agertz et al. 2006 for a discussion of numerical 
"stripping" issues). Note that while our approach does 
not distinguish between AGN formed in gas-rich "wet" 
mergers and gas-poor "dry" mergers, unlike at lower red- 
shift {e.g. Bell et aL2006), dry mergers are likely to be 
relatively unimportant at the z > 1.2 redshifts we are 
studying. 

In our fiducial run, we also assume that a fixed frac- 
tion ek = 0.05 of the bolometric energy of each AGN is 
channeled into a kinetic outflow. This value is consistent 
with other literature estimates {e.g. Furlanetto & Loeb 
2001; Nath & Roychowdhury 2002), as well as observa- 
tions (Chartas et al. 2007). Each outflow is then launched 
with an energy input of 
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Given the uncertainties surrounding AGN outflows, we 
simply model each expanding outflow as a spherical 
shell at a radius 2r v which is constructed by rearrang- 
ing all the gas within this radius, but outside r v , and 
below a density threshold of 2.5p v . Finally, the radial 
velocity v s and temperature T s of the shell are deter- 
mined by fixing the postshock temperature, T s , to be 
T s = 13.6/T [ti s /(kms -1 )] 2 and choosing v s such that the 
sum of the thermal and kinetic energy of the shell equals 
E k minus the energy used to move particles from their 
initial positions into the shell. As shown in TSC06, this 
prescription results in a level of preheating in galaxy clus- 
ters and groups that is in good agreement with observa- 
tions. 

In Figure [TJ we show the mass-averaged gas temper- 
ature in our AGN feedback and comparison runs. The 
presence of outflows has a substantial impact on the 
IGM, raising the mean gas temperature in the simulation 
by w 50%. Furthermore, as discussed in TSC06, most 
of this heating is targeted toward the densest regions, 
exactly those that contribute most to the SZ effect. 
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2.2. Construction of SZ maps 

As the thermal and kinetic SZ effect have different fre- 
quency distributions, we calculate each of them sepa- 
rately from our simulations. In the (nonrelativistic) ther- 
mal SZ case, the change in the temperature of the CMB 
as a function of frequency is given by 



AT 



y 



e x + 1 
' e x - 1 



(5) 



where x = hv/kTcMB is the dimensionless frequency, and 
the Compton y parameter is defined as 



y = a T 



dl n, 



k(T e - Tomb) 



(6) 



where ot is the Thompson cross section, m e is the mass 
of the electron, n e is the density of electrons, T e is the 
electron temperature, Tomb is again the temperature of 
the CMB, and the integral is performed over the proper 
distance along the line of sight. Note that y is simply a 
rescaled version of the line-of-sight integral of the pres- 
sure. Note also that in the Ray leigh- Jeans limit, in which 
x<l, eq. ([5]) reduces to — 2y, and we focus on this limit 
throughout our discussion below. 

For the kinetic SZ effect, on the other hand, the mag- 
nitude of the CMB distortions is frequency independent 
and AT/T = -b. In this case 



b — ctt I dln e —, 



(7) 



where v r is the radial peculiar velocity of the gas, and 
the integral is again performed along the line of sight. 

To construct maps of these distortions from our sim- 
ulations, we followed a method similar to da Silva et 
al. (2000) and Springel et al. (2002). For each simula- 
tion output, we smoothed the density, temperature, and 
velocity fields onto a cubic grid. The smoothing pro- 
cedure used the standard SPH smoothing methodology, 
that cells at a distance \v ce u — r j | from a particle i, are 
incremented by 

1 A- 

A ce u = l -miW(\v ce a - ti\,hi), (8) 

w ni Pi 

where A4 is the scalar field value for particle i, pi is the 
density of particle i, W is the B2-spline kernel, and w ni 
is a normalization factor. The smoothing process gener- 
alizes in the natural fashion for vector fields. The nor- 
malization factor is required to ensure that the weight- 
ing of the assigned kernel is correct when summed across 
a finite number of cells (the standard kernel definition 
normalizes to unity over a continuous spherical volume 
of radius 2hi). If a particle is smoothed over n ce u cells 
within a spherical volume specified by the radius 2/ij, 
then w ni is given by 



cell k 



, hi), 



(9) 



fc=i 



where we have used k to distinguish that the summation 
is over cells, rather than over particles. 

Once the smoothed grid had been constructed it was 
then projected in the x, y and z directions. Since projec- 
tion amounts to a summation along each axis direction 
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Fig. 2. — Top: Total integrated value of the evolution of the 
average Compton y parameter as a function of redshift for the 
AGN feedback run (solid line) and comparison run (dashed line) 
The shaded region is bounded below by the Scannapieco & Oh 
(2004) predictions for the y distortion due to AGN feedback (with 
efc = 0.05) and from above by the sum of this prediction plus 
the results of our comparison run. Middle: Total RMS value of 
AT/T = —2y from the thermal SZ effect, as a function of redshift, 
with lines as in the upper panel. Bottom: Total RMS AT/T from 
the kinetic SZ effect as a function of redshift. 

the smoothing process could potentially have been per- 
formed in two dimensions alone. However, we project the 
three dimensional grid to ensure that the integration in 
the projection axis accounted for the quantization effects 
associated with smoothing onto a fixed number of cells. 
Next we constructed a grid of spanning (1.1 deg) 2 and 
made up of 1024 2 rays in the fiducial run, and spanning 
(0.55 deg) 2 and made up of 512 2 rays in the comparison 
run, such that in both cases the full simulation subtended 
an angle equal to the field at our highest redshift output 
of z — 10. Finally we projected along each of our sight- 
lines, choosing a random translation and orientation at 
every redshift, and weighting the slices according to eqs. 
© and ©. 

3. RESULTS 

3.1. Overall Properties 

In Figure [H we study which redshifts contribute most 
to CMB distortions, by integrating from z = 10 (our as- 
sumed redshift of reionization) down to various redshifts. 
From eq. ((6]) the change in the average Compton y dis- 
tortion per unit dl is proportional to the mass-average 
temperature. Thus, consistent with the temperature evo- 
lution shown in Figure [TJ (y) is roughly 50% higher in 
the AGN feedback run at all redshifts. Furthermore, for 
most redshifts, this increase is similar to the one derived 
from analytic estimates as presented in Scannapieco &z 
Oh (2004). In fact, the difference between this model and 
our simulations may be due to the presence of small-scale 
structure around AGN outflows, which were ignored in 
the analytic estimates, and whose mixing with outflow 
material may be somewhat under-resolved in the simu- 
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Fig. 3. — Map of the thermal SZ effect in our AGN feedback 
and comparison runs (large panel and inset panel, respectively). 
Here the area of the largest panel is (l.ldeg) 2 and the color scale 
runs from ^£ = (black) to -5 X 10~ 5 (white), (or to -135 
IxK). 

lations (see TSC06 for a full discussion of this issue). 

Yet, the average SZ distortion is not the most eas- 
ily detectable quantity. Rather, the majority of CMB 
measurements are differential in nature, and sensitive to 
small changes on top of an overall background that is 
much less well measured. To quantify these differences, 
in the center and low panels of Figure 2, we plot the 
RMS scatter in AT/T due to the thermal and kinetic 
SZ effects respectively. In both cases the changes in the 
amplitude of the distortions are small, indicating that 
AGN have a minor impact on the overall variance of the 
temperature and velocity. 

In Figure [3] we directly compare maps of the thermal 
SZ distortions from both these runs, integrated down to 
our final simulation rcdshift of z = 1.2. The maps clearly 
indicate a lack of small hot regions in the AGN feedback 
run relative to the comparison run. To help bring out the 
structure in these maps, we also plot the same data on a 
logarithmic scale in Figured! Here we see that the AGN 
feedback run has a higher average level of distortions, 
consistent with the (y) evolution in Figure [2j Further- 
more, while it is difficult to tell the overall level of the 
variance between the maps, it is clear that the scale of 
the structures is quite different. In particular, while the 
AGN feedback simulation has fewer small pockets of hot 
gas, these are compensated for by a number of larger and 
more diffuse heated regions. 

Finally, the kinetic SZ maps, as shown in Figure [5) 
look similar in the two simulations. The differences 
in this case appear to be quite subtle, and must be 
distinguished by a more detailed, quantitative approach. 
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Fig. 4. — Map of the logarithm of the thermal SZ ef- 
fect, ranging from log(-AT/T) = log(l X 10~ 6 ) (black) to 
log(— AT/T) = log(5 X 10 -5 ) (white), with panels as in Figure 
[3] While the AGN simulation has a higher average Compton 
distortion, structure has been smoothed out on the smallest scales. 




Fig. 5. — Map of the kinetic SZ effect ranging from 
■^f = —1 X 10~ 5 (black, moving away from the observer) to 
1 X 10~ 5 (white, moving towards the observer), that is -27.3 to 
27.3 /uK. 
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3.2. Power Spectrum 
3.2.1. Results from Simulated Maps 

The most common measure of CMB fluctuations is 
the angular power spectrum, obtained by decomposing 
the temperature distribution on the sky into spherical 
harmonics, AT/T(n) = Y^i, m a t,mYt,m(&-), and carry- 
ing out an average over the coefficients, Ci = (2£ + 
1) _1 J2m=-t. a i,mO*t m - For a small field of view, as is the 
case for our constructed map, this is equivalent to per- 
forming a Fourier Transform to obtain 



AT 
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d 2 9 exp 



(-a.e)^L(6), (io) 



(where and k are two dimensional vectors in the plane 
of the sky) and then carrying out an azimuthal average 
to obtain 



C t = P a „g(K = t) = 7^ 
Ztt 



AT* 



{k,4>), 



where k — (ncostp, /tsin<^>). Finally, as an additional 
check one can use the fact that the mean squared tem- 
perature variance is equal to ^2 e (2£ + 1)0^/(4%). 

In Figure [6] we plot the angular power spectrum from 
our AGN feedback and comparison simulations, averag- 
ing over 16 realizations of the maps: (4.4 deg) 2 in the 
AGN feedback case and (2.2 deg) 2 in the comparison 
case. As expected from a visual inspection of Figures [3] 
and S], AGN feedback has a clear impact on the thermal 
SZ signal, smoothing the structures at the smallest scales 
(£> 10 4 ) while at the same time adding large heated re- 
gions that increase the overall fluctuations at larger scales 
(£ < 10 3 ). On the other hand, AGN feedback has only 
a weak impact on the kinetic SZ effect, slightly smooth- 
ing the smallest structures, without having a noticeable 
impact on larger scales. To some extent, the lack of this 
large-scale signal may be due to our choice of initially 
spherical outflows, which should not significantly add to 
the total integrated column of material moving towards 
or away from the observer along any line of sight. If out- 
flows were extremely asymmetric, however, from eqs. ([6]) 
and ([7]) the ratio of the distortions in the kinetic to the 
thermal SZ effects would go as 



v/c 



y 



kT/m e c 2 m p (v/c) 



(12) 



For a typical outflow velocity of v 1000 km/s, this 
gives b/y ~ 1/6, such that the even in this extreme case, 
the increase in Cis due to the kinetic effect would be at 
least ss 30 times smaller than the thermal signal. 

However, even in the thermal case, the increase in the 
large-scale angular power spectrum is not dramatic. For 
comparison, in Figure[5]we have also plotted the primary 
CMB signal, as well as a summary of recent measure- 
ments from the CBI and ACBAR experiments. From 
these points it is clear that the secondary anisotropics 
seen in current experiments are not likely to have been 
imprinted by SZ distortions at z > 1.2, nor does the pos- 
sible excess of small-scale power seen in the CBI experi- 
ment (see however, Kuo et al. 2007) appear to be related 
to AGN feedback. 
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Fig. 6. — Top: Power spectrum of the thermal SZ distortions 
in the AGN feedback run (thick solid line) and comparison run 
(thick dashed line), integrated down to Z = 1.2. The open squares 
are measurements by the CBI experiment (Bond et al. 2005), the 
solid points are measurements by the ACBAR experiment (Kuo et 
al. 2007), and the stars are measurements from the BIMA CMB 
Anisotropy Survey (Dawson et al. 2006). These are plotted around 
the primary anisotropy signal (dotted lines). For comparison, the 
thin solid and dashed lines show the kinetic SZ anisotropies for 
the AGN feedback and comparison run, as in the bottom panel. 
The dot-dashed lines show the level of anisotropies detectable 
by an experiment spanning a (20 deg) 2 area with a $fwhm = 2 
arcmin beam, with a noise level of (10 /iK) 2 per arcmin 2 (upper 
line) and with a 6>p-yvHM = 1 arcmin with a noise level of (2 
/aK) 2 per arcmin 2 (lower line). Bottom: Power spectrum from 
the kinetic SZ distortions in the AGN feedback run (thick solid 
line) and comparison run (thick dashed line), integrated down to 
z = 1.2. Points and dot-dashed lines are as in the upper panel, and 
the thin lines are the thermal SZ distortions, shown for comparison. 

3.2.2. Required Sensitivity and Resolution to detect AGN 

distortions 

Pushing to the near future, however, the increase in 
the angular power spectrum due to AGN feedback is well 
within the range of upcoming experiments. To calculate 
these sensitivities we make use the simple estimate from 
Jungman et al. (1996), which gives the RMS noise at mul- 
tipole i as 



1-1/2 



(2*+l)/*y 



Ci + (W/ S ky)~ 



(13) 



where w = (ffpix^FWHM) -2 is the weight per solid an- 
gle, a pixel-size independent measure of the noise that is 
calculated from the full-width at half-max of the Gaus- 
sian beam, #fwhm, and the noise variance per #fwhm x 
^fmwh pixel, CT 2 ix . Finally, a h = 7.42 x 10 -3 x # F hwm 
is a measure of the beam size and / s k y is the fraction of 
the sky covered. 

Comparing these measurements with the Ces from our 
simulations, we see that the excess due to AGN feedback 
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is easily measurable by an experiment with a typical full- 
width at half max beam size of #fwhm ~ 1 arcmin and a 
noise level of about (2 fJ,K) 2 per arcmin 2 , scanning over a 
(20 deg) 2 patch of the sky. However, the most convincing 
constraints on AGN feedback are not likely to come from 
statistical constraints derived purely from the microwave 
background. 

3.3. Cross Correlations 

Unlike the primary microwave background signal, the 
SZ distributions are not well described by a Gaussian 
random field. Rather they contain a wealth of informa- 
tion beyond the angular power spectrum. In fact, as is 
particularly clear from Figures [3] and [H the thermal SZ 
distribution is much more reminiscent of lower-redshift 
galaxy surveys than z « 1000 measurements. Further- 
more, the structures seen in these maps are imprinted 
by galaxy activity, meaning that they are strongly corre- 
lated with the positions of sources in large-field optical 
and infrared surveys. With this in mind, we considered 
the cross-correlation of our maps with three types of op- 
tical sources: 

• quasars, identified as newly-formed black holes that 
shine for t& — 0.055r v /v c of the dynamical time of 
the system (as in §2.1) 

• post-starburst elliptical (E+A) galaxies, identified 
as mergers observed within 200 Myrs of coalescence 

• elliptical galaxies, identified as mergers observed at 
any time after coalescence. 

In Figure [7] we show the results of coadding the ther- 
mal SZ distortions around each of these objects as se- 
lected in a (2.2 deg) 2 region, made up of 4 maps in the 
AGN feedback run and 16 maps in the comparison run. 
The most obvious feature in this plot is the overall higher 
level of the SZ background in the AGN feedback run, a 
feature that is difficult to observe directly, as discussed 
in §2.1. Beyond this overall offset however, these images 
also uncover a substantial "halo" of SZ distortions, which 
is much higher in amplitude and more spatially extended 
in the AGN feedback case. While this excess is somewhat 
lost in the noise for quasars, which are the rarest sources, 
it is clearly and dramatically present in the coadded maps 
of the much more numerous E+A and elliptical galaxies. 
As expected from our discussion above, constructing sim- 
ilar coadded images of the kSZ effect yielded differences 
that were at least an order of magnitude smaller. 

To further study the detectability of the much more 
promising thermal signal we processed each of the im- 
ages in Figure [7] in a manner similar to how one might 
work with real data sets in the near future. First, we 
shifted each of them by the overall average y distortion 
to account for the difficult to observe offset between the 
AGN feedback and comparison cases. Next, to approx- 
imate the angular resolution of the next generation of 
CMB telescopes, we convolved each image with Gaus- 
sian beams with #fwhm of 1 and 2 arcmins. Finally, as 
coadding sources washes out any rotational asymmetries, 
we performed an azimuthal average to reduce noise while 
retaining the same information. 

The results of this procedure are shown in Figure [5J 
for quasars with B-band luminosities above 1 x 10 13 Lq, 



3 x 10 12 L©, and 1 x 10 12 L Q , which correspond to black 
hole masses of 3 x 10 9 M Q , 1 x 10 9 M Q , and 3 x 1O 8 M . 
Each panel in this plot shows not only the average signal 
over all sources in a (2.2 deg) 2 region of the sky, but the 
maximum and minimum signal in any of the four (1.1 
deg) 2 quadrants in this region, to give an impression 
of the variation in this quantity across the sky. While 
signals can vary substantially between quadrants, in all 
cases the overall average SZ signal is significantly higher 
in the AGN feedback case. As an estimate of the sensi- 
tivity of future experiments we also include on this plot 
a background noise level calculated as 



N. 



sources 



(2n9d6) 



1/2 



(14) 



where d9 is the 1.1 deg/1024 = 0.0644 arcmin pixel size 
of our maps, N SOUIces is the number of sources that would 
be found in a (20 deg) 2 region of sky, and the pixel noise 
level, (7pi x , is again estimated as (10 /^K) 2 per arcmin 2 
(cross-hatched region) and (2 /xK) 2 per arcmin 2 (solid 
region) . 

The impact of AGN outflows is well above the noise for 
all sources, even if 6>fwhm = 2' and cr 2 ix = (10 ^K) 2 per 
arcmin 2 . However, in most cases, the profile of the SZ 
signal above the noise is indistinguishable from a point 
source (which would appear on this plot as an inverted 
parabola that drops by a factor of two at a distance of 
#fwhm/2). At the same time quasars themselves are 
often significant microwave point sources. To estimate 
this intrinsic contribution we convert the flux per unit 
frequency, F v , to CMB temperature units as 



f dB v 
V dT 



9 2 

FWHM 



0.0084 mJy \ 6» F whm / 



1 arcmin \ (e x — 1) 



F v - 



where B v is the Planck Function and, as in eq. (|SJ), 
x = hv/kTcMB (Scott & White 1999). The ratio of 
vL v at microwave wavelengths to the quasars bolometric 
luminosity is ~ 1/1000 for "radio loud" objects which 
make up about 10% of the population, and several or- 
ders of magnitude less for radio-quiet objects (Elvis et 
at 1994). Estimating a typical luminosity distance to a 
quasar as « 3000 Mpc and considering a typical ob- 
serving frequency of 100 GHz this gives 



AT: 



Mbh 
10 6 M P 



1 arcmin \ ^ 

#FWHM / 



(16) 



which is much larger than the thermal SZ signal. 

While this can be reduced substantially by selecting 
radio-quiet quasars, a better approach is to work with 
objects observed at later times after the merger. In Fig- 
ure [H] we show the azimuthally averaged excess thermal 
SZ contribution around post-starburst galaxies, identi- 
fied as mergers observed after the quasar phase, but 
within 200 Myrs after coalescence. Here we have con- 
verted each black hole mass to its corresponding stellar 
bulge mass using a factor of w 400 as obtained from 
the analysis in Marconi & Hunt (2003). Working with 
these sources has two main advantages. Firstly, as they 
are much more common, one can consider slightly larger 
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Fig. 7. — Postage stamp images of the thermal SZ effect in a 6 X 6 arcmin 2 region around various objects. Each pair of panels 
shows the distortions in the AGN run (left) and comparison run (right) averaged over all the sources in a (2.2 deg) 2 region of the sky. 
From top to bottom, the rows are labeled by the black hole masses, which range from 1O 1O M0 to 3 X 1O 8 M0. From left to right each 
of the pairs correspond to quasars, post-starburst ellipticals, and all elliptical galaxies. All panels are shown using a log scale from 
log(-AT/T) = log(2 x 10" 5 ) (white) to log(-AT/T) = log(2 x 10" 6 ) (black). 



bulges and hence larger and more spatially extended SZ 
distortions, while at the same time improving the num- 
ber of sources that can be coadded in a fixed region of the 
sky. Furthermore, as the nuclear luminosities are typi- 
cally more than 100 times less during this phase (e.g. 
Hopkins et al. 2007) , point source subtraction is much 
easier, although some exclusion of radio-loud sources may 
still be necessary. 

In Figure [TOl we show the SZ contribution obtained by 
coadding all bulges observed at any time after the initial 
quasar phase. Again we have converted black hole mass 
into bulge mass as well as into overall B-band luminosity 
following observed relations (Marconi & Hunt 2003). Ex- 
tending our analysis in this way has only a weak impact 
on the SZ profile, while at the same time improving the 
statistics even further. In this case the complete spatial 
profiles of the sources are observable with more than 10 
a precision out to very large radii, clearly tracing out 
the impact of AGN feedback. Note that these sources 
are relatively easy to detect, as even given a typical lu- 
minosity distance of 10 4 Mpc, they are readily observed 
by a photometric survey with an overall magnitude limit 
ofm s «22. 

Furthermore, the level of thermal y distortions ob- 
served by coadding around bulges is a precise and linear 
measure of the most important quantity for modeling the 
impact of AGN outflows on galaxy formation. Integrat- 
ing eq. © over a patch of sky around a source we have 



d9y{9) 



Ti. e C 2 ^ang 



dVn e (V)k[T e (V)-T CM B], 

(17) 



where l ang is the angular diameter distance to the source, 

9 is a vector in the plane of the sky in units of radians, 
and the volume integral is performed over the heated 
region around the source. But this integral is simply 
(2/3)£thermai(l + A) / (2 + A) where A = 0.08 is the cos- 
mological number abundance of helium. The linc-of-sight 
integral of the pressure has thus been transformed into 
a volume integral of the pressure. This means that by 
integrating the Compton distortions over the sky we can 
directly measure the thermal energy added to the IGM 
for each object. Rewriting eq. Q17[) in terms of angles in 
arcminutes and the low-frequency microwave background 
distortions in units of \xK this becomes 



thermal 



-4.8 x 10"° ergs p f? AT % (18) 



fiK arcmin 



where l ang is the angular distance in units of 3000 Mpc 
and -Ethcrmai is the total excess thermal energy associated 
with the source, that is the thermal energy gained from 
the initial collapse of the baryons, plus the contribution 
from the AGN, minus losses due to cooling and PdV 
work done during expansion. 

3.4. Measuring AGN Feedback 

Putting these results together, a clear and direct way 
to constrain AGN feedback presents itself. This can be 
summarized as follows: 

• Obtain SZ data with 1 or 2 arcminutc resolution 
at a noise level of a few fiK per arcmin 2 in a ~ 
400 deg 2 patch of the sky in which near infrared 
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Fig. 8. — Average excess thermal SZ distortions around quasars 
in our simulations. In each panel the thick solid line shows AT 
in the AGN feedback run, averaged over a (2.2 deg) 2 area of 
the sky. As a measure of the variance over the sky, the thin 
solid line shows the maximum and minimum AT signal in single 
(1.1 deg) 2 quadrants of this (2.2 deg) 2 map. Similarly, the thick 
dashed line shows the average AT in (2.2 deg) 2 of sky from 
the comparison run, bracketed by the maximum and minimum 
signals in quadrants. Finally, the shaded region shows the noise 
level attainable by averaging over all sources in (20 deg) 2 , for 
an experiment with a noise-per-pixel of (10 /iK) 2 per arcmin 2 
(cross-hatched regions) and (2 /iK) 2 per arcmin 2 (solid regions). 
The left panels show distortions as observed with a Sfwhm = 1' 
beam and the right panels show the same distortions for a 
6*FWHM = 2' beam. From top to bottom the panels correspond 
to black hole masses of 3 X 10 9 A4q, 1 X 10 9 M Q , and 3 X 10 s M®, 
respectively, as marked on the left panels, and the equivalent 
B-band luminosity is given in the right panels. 

or optical photometry is available to a moderately 
deep magnitude limit (such as nis ~ 22). 

• Select all quiescent early-type galaxies with pho- 
tometric redshifts in the essential z = 1 — 3 red- 
shift range for AGN feedback. Alternatively, if suf- 
ficiently good photometry is available, one can also 
study the subset of post-starburst (E+A) galaxies. 

• If possible, use radio observations to reject point 
sources that may contaminate the thermal SZ sig- 
nal. 

• Compute the excess thermal SZ signal summed 
within a sa 2 arcminute radius around each source, 
and use this to compute the total thermal energy 
as per eq. (fT8|) . 

• Bin the results as a function stellar mass, and av- 
erage over each bin to compute the total average 
IGM thermal energy input associated with bulges 
as function of their stellar mass. 

The results of carrying out this procedure over (2.2 
deg) 2 of sky calculated from our simulations are shown 




arcmin 

Fig. 9. — Average excess thermal SZ distortions around 
post-starburst galaxies in our simulations, observed within 200 
Mys of a merger. As there are many more of these objects than 
quasars, we select slightly larger objects in this figure than in 
Figure [8] Thus from top to bottom the panels correspond to black 
hole masses of 1 X 10 10 M Q , 3 X 10 9 Af Q , and 1 X KFM Q (marked 
on the left panels), or equivalent bulge masses of 4 x 10 12 M Q , 
1 X 10 12 M Q , and 4 X 10 n Af Q (marked on the right panels). 
Otherwise the panels, lines, and shaded regions are defined as in 
Figure [8] Thus, the left panels show distortions as observed with a 
$FWHM = 1' beam and the right panels show the same distortions 
for a Sfwhm = 2' beam; The thick solid (AGN feedback) and 
dashed (comparison) lines show the average value of the SZ 
distortions averaged over all sources in a (2.2 deg) 2 region of the 
sky, bracketed by the maximum and minimum signal in (1.1 deg) 2 
quadrants; and the shaded regions show the noise levels attainable 
by averaging over (20 deg) 2 , with instrument noise as in Figure [8] 

in Figure 1111 Here we have grouped all bulges in loga- 
rithmic bins of stellar mass with width 0.25, again taking 
a fixed ratio of 400 between Mb u i ge and Mbh- To cal- 
culate the uncertainty in each bin we have accounted 
both for Poisson_ noise and intrinsic scatter, taking 

f^bin (-^thermal) — [-^thermal, bin ~t~ ^"source (-^thermal)]/ V^Vb in ? 

where -Ethermai.bin is the average thermal energy in a 
given bin, <r S ource (-^thermal) is the RMS scatter between 
sources in a given bin, where we computed A^in conser- 
vatively as the number of sources in a given bin relative 
to a signal realization of the map (1/4 of the [2.2 deg] 2 
region in the AGN feedback run and 1/16 of this region 
in the comparison run). 

At both high and low redshifts the AGN feedback run 
shows a clear excess, which scales linearly with Mbuige 
as expected from eq. (jj). For comparison, we use this 
equation to plot the energy added to the IGM as a func- 
tion of bulge mass. Since kinetic energy is converted into 
thermal energy as the outflows accrete material, and as 
radiative losses are small for these objects (e.g. Oh & 
Benson 2003; Scannapieco & Oh 2004), we expect that 
most of this energy should be observable as -Ethermai ■ In 
fact, this figure shows that at both low and high redshift, 
the thermal SZ excess closely traces this energy input as 
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Fig. 10. — Average excess thermal SZ distortions around all 
bulges in our simulations. Mass limits, beam sizes, and symbols 
are as in Figure [9] Right panels are labeled with estimated total 
bulge Lb luminosities following Marconi & Hunt (2003) 
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Fig. 11. — IGM thermal energy computed from eq. I|18jl as a 
function of bulge stellar mass in our AGN feedback (solid points) 
and comparison (crosses) simulations. In each bin the point shows 
the results of averaging over all sources in a logjQ(Mb u i ge ) = 0.25 
bin surrounded by error bars that account for both the intrinsic 
variation between sources and the fact that multiple maps have 
been derived from the same simulation (see text). The solid lines, 
on the other hand, show the energy input added around bulges in 
our AGN feedback simulation, as per eq. JJ). In the upper panel 
we consider all bulges in the redshift range from z = 1.2 — 2.0, and 
in the lower panel from z = 2.0 — 3.0. 
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Fig. 12. — Thermal energy around post-starburst galaxies as 
a function of bulge stellar mass. Panels and symbols are as in 
Figure [TTJ 

a function of bulge mass. Thus AGN feedback is not only 
detectable by the SZ effect, but the level of this feedback 
as a function of mass can be obtained directly from SZ 
measurements . 

Finally, in Figure [T^l we plot thermal for the post- 
starburst galaxies in our simulations. As there are fewer 
of these objects, this results in a higher overall scatter 
than in the case of all bulges, but the results are other- 
wise similar and consistent with the level of excess en- 
ergy added to the AGN feedback run. While this is to 
be expected given the instantaneous and impulsive na- 
ture of the outflows in our simulation, comparisons of 
post-starburst with other bulges may nevertheless help 
to distinguish this type of feedback from more gradual 
AGN heating, as suggested, for example, in Croton et 
al. (2006). 

4. CONCLUSIONS 

Microwave background measurements have played a 
key role throughout cosmology: establishing the mod- 
ern inflationary model, constraining the initial conditions 
for structure formation, and providing precise measure- 
ments of cosmological parameters. Yet at smaller scales, 
the potential of CMB measurements is largely untapped. 

At the same time, detailed observations at lower red- 
shifts have uncovered the central importance of nongrav- 
itational heating, which is likely to be associated with 
outflows from AGN. Although the details of this pro- 
cess remain uncertain, we have shown using a suite of 
large numerical simulations that the key observations 
constraining this process may once again come from the 
microwave background. While AGN outflows impact the 
CMB both through their peculiar motions and through 
IGM heating, it is the thermal SZ signal that provides 
the most detailed measurements. Working purely with 
CMB observations, AGN feedback can be constrained 
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both through a difficult to measure overall offset in the 
frequency spectrum, and a more easily detected shift in 
CMB anisotropics from rs 1/10 arcmin to w 1 arcmin 
scales. 

Although changes in the angular power spectrum from 
this shift are marginally detectable with the upcom- 
ing generation of small-scale microwave telescopes, the 
true test of AGN feedback will come from the cross- 
correlation of this data with infrared and optical surveys. 
In particular, coadding the CMB signal around quasars, 
post star burst, and quiescent elliptical galaxies from our 
simulations shows a systematic "halo" of SZ distortions 
around each of these objects, which is much higher in 
amplitude and more spatially extended in the AGN feed- 
back case. Also, while contamination from point source 
emission make quasars difficult to work with, the excess 
signal is clearly detectable for both E+A and older el- 
liptical galaxies, which can be selected with only moder- 
ately deep tub ~ 22 photometry. 

Furthermore, as the thermal SZ effect is proportional 
to the line of sight integral of the pressure, summing up 
the excess distortion in a patch of sky around each type 
of source provides a sensitive and direct measure of the 
thermal energy associated with feedback. In fact, carry- 
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ing out this procedure on our simulation we have shown 
that we can easily recover not only the overall level of 
AGN feedback, but its dependence on galaxy mass, red- 
shift, and galaxy type. Although it provides a popu- 
lar and elegant solution to many outstanding problems, 
the impact of IGM heating by AGN outflows remains 
perhaps the most uncertain outstanding issue in galaxy 
formation. Working together with large-field optical sur- 
veys, small-scale CMB experiments will soon be able to 
place strong constraints on this missing piece in our phys- 
ical understanding of the history of our universe. 
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